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ABSTRACT 

We study the evolution of the mass autocorrelation function by describing 
the growth of density fluctuations through the Zel'dovich approximation. The 
results are directly compared with the predictions of the scaling hypothesis 
for clustering evolution extracted from numerical simulations (Hamilton et 
al. 1991), as implemented by Jain, Mo & White (1995). We find very good 
agreement between the correlations on mildly non-linear scales and on com- 
pletely linear scales. In between these regimes, we note that the density fields 
evolved through the Zel'dovich approximation show more non-linear features 
than predicted by the scaling ansatz which is, however, forced to match the 
linear evolution on scales larger than the simulation box. In any case, the 
scaling ansatz by Baugh & Gaztahaga (1996), calibrated against large box 
simulations agrees better with ZA predictions on large scales, keeping good 
accuracy also on intermediate scales. 

We show that mode-coupling is able to move the first zero crossing of 
£(r) as time goes on. A detailed fit of the time dependence of this shifting is 
given for a CDM model. The evolution of the cross correlation of the density 
fluctuation field evaluated at two different times is also studied. The possible 
implications of the results for the analysis of the observed correlation function 
of high redshift galaxies are discussed. 

Key words: galaxies: clustering - cosmology: theory - large-scale structure 
of Universe 

1 INTRODUCTION 

In the last two decades, redshift surveys provided a wealth of informations about the spatial distribution of local 
galaxies, revealing the existence of large-scale structures. The most widely used statistical tool to quantify the degree 
of clustering has been the galaxy two-point correlation function, both in its angular w g (6) and spatial £ 9 (r) versions 
(see, e.g., Peebles 1980). As previous studies were confined to the nearby universe, nowadays new observational 
resources permit to extend the correlation analysis to deeper samples. In fact, the Canada-France Redshift Survey 
has recently provided the new opportunity to investigate the clustering properties of galaxies out to redshifts z ~ 
1 (Le Fevre et al. 1996). Moreover, only lately it has been possible to analyse the angular distribution of faint 
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galaxies by using the Hawaii Keck K-band survey (Cowie et al. 1996) and the Hubble Deep Field data (Villumsen, 
Freudling & da Costa 1997). Therefore, we are now confident of detecting a direct signature of redshift dependence 
in the observed correlation function. For this reason, the theoretical analysis of the evolution of the mass two-point 
correlation function, £(r), is becoming a fundamental topic of modern cosmology. However, it is worth stressing that 
the interpretation of the observational data is not immediate: before obtaining the "real" change of the large-scale 
structure one has to consider the possible evolution of the galaxy population (as well as the related selection effects) 
and of the bias factor that formally relates £ 9 to £ (see, e.g., Matarrese et al. 1997). 

The observational results should then be compared to the predictions of the existing models for structure for- 
mation. One of the several issues involved in this comparison is represented by the lack of a standard description of 
clustering evolution: analytical treatments are generally unable to manage this fully non-linear problem while numer- 
ical simulations are limited in resolution. However, new light has been recently shed on this argument. Hamilton et 
al. (1991) suggested that the correlation function obtained through N-body simulations of an Einstein-de Sitter uni- 
verse, in which the structure develops hierarchically, can be easily reproduced by applying a non-local and non-linear 
transformation to the linear £(r). This ansatz has been refined and extended to more general cosmological scenarios 
by a number of authors (Peacock & Dodds 1994, Jain, Mo & White 1995, Peacock & Dodds 1996). Moreover it is 
possible to give theoretical arguments that account for the scaling hypothesis (Nityananda & Padmanabhan 1994). 

The main purpose of this paper is to compare the predictions of the Zel'dovich approximation (Zel'dovich 1970) 
with the scaling ansatz formulated in the version of Jain, Mo & White (1995, hereafter JMW). Actually, it would be 
very interesting to obtain all the details of the semi-empirical scaling relationship in the framework of the gravitational 
instability scenario. However, in the absence of a model for the advanced phases of clustering evolution, we are forced 
to analyse only the onset of non-linear dynamics. 

From the theoretical point of view, the evolution of the two-point correlation function is strictly related to 
the dynamical development of the density field g(x,i). When the dimensionless density contrast <5(x, t) = [(?(x,i) — 
£(*)]/(?(*) is much smaller than unity, the growth of the fluctuations can be followed performing a perturbative 
approach (see, e.g., Peebles 1980, Fry 1984, Scoccimarro & Frieman 1996a). At the lowest order (linear theory), the 
different Fourier modes of <5(x,t) evolve independently provided that their power spectrum is less steep than k 4 at 
small k (Zel'dovich 1965, Peebles 1974). As the fluctuations grow, however, the interactions between different modes 
become more and more important. The effect of this mode-coupling on the two-point statistics has been studied 
by many authors using higher-order than linear terms in the perturbation expansion. Juskiewicz, Sonoda & Barrow 
(1984) computed the second-order contribution to £(r) for an exponentially smoothed linear spectrum P(k) tx k 2 , 
finding that non-linear interactions among long wavelength modes act as a source for short A perturbations. As a 
matter of fact, they found a substantial decrease of the characteristic scale of clustering with the evolution. However 
Suto & Sasaki (1991) and Makino, Sasaki & Suto (1992), analysing exponentially filtered scale-free spectra, found 
that second-order effects can either suppress or enhance the growth of perturbations on large scales, depending on the 
shape and the amplitude of the fluctuation spectrum. In particular Makino, Sasaki & Suto (1992), modelling a CDM 
spectrum with two different power laws, concluded that the effects of mode-coupling are generally very small and 
completely negligible on scales rjj 20ft~ 1 Mpc (where h denotes the Hubble constant in units of 100 kms -1 Mpc -1 ). 
The second-order correction to the "true" linear CDM spectrum has been calculated by Coles (1990) who computed 
also the respective correlation function. The results show that, for moderate evolution, the large-scale distortions are 
of no importance, while later (for ct 8 ~ 1, where as represents the rms linear mass fluctuation in spheres of radius 
8/i _1 Mpc) non-linear effects can increase the clustering strength on scales r > 35 ft x Mpc; for example, the first 
zero crossing of £(r) can be significantly shifted with respect to linear predictions. Similar results were obtained by 
Baugh & Efstathiou (1994) who also found good agreement with the output of numerical simulations. However, Jain 
& Bertschinger (1994) pointed out that the perturbative approach is able to reproduce the N-body outcomes only 
at early times (agj; 0.5 — 1). Moreover, the recent analysis applied to scale-free spectra by Scoccimarro & Friemann 
(1996b) showed that the validity of perturbation theory is restricted to a small range of spectral indices. 

In this paper, we want to study the non-linear evolution of the mass autocorrelation function by describing the 
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growth of density fluctuations through the Zel'dovich approximation (hereafter ZA). In effect, Eulerian second-order 
perturbation theory may break down once the mass variance becomes sufficiently large. On the other hand, we know 
that ZA, especially in its "truncated" form, is able to reproduce fairly well the outcomes of N-body simulations even 
in the mildly non-linear regime (Melott, Pellman & Shandarin 1994). The main advantage of ZA over other dynamical 
approximations (for a recent review see, e.g., Sahni & Coles 1995) is that it permits analytical investigations ensuring 
at the same time good accuracy, at least for quasi-linear scales. The pioneering analysis by Bond & Couchmann 
(1988) showed that ZA is able to predict the shifting of the first zero crossing of the correlation function. In Section 3 
we will give a detailed quantitative description of this effect. Other features of the mass two-point correlation function 
in ZA have been discussed by Mann, Heavens & Peacock (1993, hereafter MHP). Moreover, the related evolution 
of the power spectrum has been studied by Taylor (1993), Schneider & Bartelmann (1995) and Taylor & Hamilton 
(1996). These authors showed that ZA is able to describe the generation of small-scale power through mode coupling, 
at least at early times. Besides Fisher & Nusser (1996) and Taylor & Hamilton (1996) succeeded in computing the 
power spectrum also in redshift space. 

This paper is organized as follows. In Section 2 we briefly introduce the Zel'dovich approximation while in Section 
3 we compute the cross correlation function between the mass density field evaluated at two different times. The usual 
two-point correlation function is obtained as a particular case of this more general quantity. The redshift evolution of 
£(r) in a CDM model is the last subject of Section 3. In Section 4 we compare the predictions of ZA with the scaling 
ansatz of JMW. In Section 5 we use our results to evaluate the correlation function of a collection of objects sampled 
by an observer in a wide redshift interval of his past light cone. We then propose a simplified scheme to compute this 
quantity so as to improve another approximation presented in the literature. A brief summary is given in Section 6. 



2 THE ZEL'DOVICH APPROXIMATION 

Let us consider a set of collisionless, self-gravitating particles in an expanding universe with scale factor a(t). We can 
describe the motion of each point-like particle writing its actual (Eulerian) comoving position, x, at time t as the 
sum of its initial (Lagrangian) comoving position, q, plus a displacement: 

x(q,t) = q+S(q,t). (1) 

The displacement vector field S(q, i) represents the effect of density perturbations on the trajectories. The Zel'dovich 
approximation is obtained by assuming the separability of the temporal and spatial parts of S(q, i) and by requiring 
equation ([!]) to give the correct evolution of <5(x, t) in the linear regime. Considering only the growing mode for a 
pressureless fluid, one gets (Zel'dovich 1970): 

S(q,i) = -b(t)V0| q (2) 

where b(t) is the linear growth factor and <^>(q) represents the initial peculiar velocity potential that at the linear 
stage is proportional to the gravitational potential <I>o(q). The Zel'dovich approximation can be also extracted from 
a fully Lagrangian approach to the evolution of density fluctuations (Buchert 1989, Moutarde et al. 1991, Bouchet 
et al. 1992, Buchert 1993, Catelan 1995). In this case, ZA corresponds to the first order solution provided that the 
initial velocity field is irrotational and the initial peculiar velocity and acceleration fields are everywhere parallel. 

Equations (|l]) and (Q) define a mapping from Lagrangian to Eulerian space that develops caustics as time goes 
on (Shandarin & Zel'dovich 1989). However, the "Zel'dovich fluid" is a system with infinite memory: even after the 
intersection of two trajectories, the motion of the particles is determined by their initial conditions according to 
equation (^) . The lack of self-gravity between intersecting streams causes the forming structure to be rapidly washed 
out. This is a severe problem especially in hierarchical models of structure formation, where caustics appear early on 
small scales causing ZA to become soon inaccurate. Nevertheless Coles, Melott & Shandarin (1993) showed that a 
modified version of ZA, the "truncated" ZA, obtained by smoothing the initial conditions, is able to reproduce with 
good accuracy the density distributions obtained from numerical simulations. Melott, Pellman & Shandarin (1994) 
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found that the optimal version of the truncation procedure is accomplished by using a Gaussian window to smooth 
the linearly extrapolated power spectrum of the density fluctuation field b 2 (t)P(k): |J] 

Pr(k, t) = b 2 (t)P(k) exp [-k 2 R 2 f (t)] (3) 

where the filtering radius Rf(t) increases with time being related to the typical scale going non-linear. The success 
of this approximation can be justified by noticing that the non-linearly evolved gravitational potential resembles its 
smoothed linear counterpart (Pauls & Melott 1995). In the following we will adopt the filtering prescription given in 
equation (J3f) - 



3 THE TWO-POINT CORRELATION FUNCTION IN THE ZEL'DOVICH APPROXIMATION 



Assuming that initially the mass is evenly distributed in Lagrangian space, implies that the Eulerian density field is 
related to the Lagrangian displacement field via the relation: 



g(x,t) = e(t) / d 3 q5 D [x-q-S(q,t)], 



(4) 

where 5jj(x) denotes the three-dimensional Dirac delta function. For purposes that will be clarified in Section 5, we 
are interested in computing the cross correlation function between the density contrast field evaluated at two different 
times: 

<<5(xi,ii)<5(x2,i 2 )} = { J d 3 qi d 3 q 2 8 D [xi -qi -S(qi,ii)]fo [x 2 - q 2 - S(q 2 , t 2 )]) - 1 (5) 

where {•) represents the average over an ensemble of realizations. Before going any further, it is convenient to Fourier 
transform the Dirac delta functions in equation (0) obtaining: 



1 + (<5(xi,ti)5(x 2 ,t 2 )) = / tfgitfga 



,;; d 3 Wl d 3 W 2 



(2tt) 3 (2tt) : 



■ exp 
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(6) 



We then use equation (J2p to introduce ZA into equation (g). In such a way, by assuming, as usual, that <^>(q) is a 
statistically homogeneous and isotropic Gaussian field, uniquely specified by its power spectrum P<f,(k) oc P(k)/k 4 , 
the ensemble average contained in equation (^) can be written as a functional integral: 
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(7) 



where the kernel A"(q, q') represents the functional inverse of the two-point correlation function of the field 0(q). By 
defining a six-dimensional vector c* = (wi, w 2 ) and choosing the z-axis of our reference frame in the direction of the 
vector q = qi — q 2 , we can reduce equation ([?]) to the form: 



(exp 



w £ ■ S(q/,t/) 



> = exp 
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(«) 



where the matrix M has the structure 
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t We set b = 1 at the present epoch. 
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Figure 1. Left panel: the mass autocorrelation function, obtained using ZA, for a COBE normalized CDM linear spectrum 
is plotted for different values of the truncation radius Rf (in /i _1 Mpc). Right panel: dependence of the correlation function 
evaluated at r = lh~ 1 Mpc on Rf. 
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having denoted by ji(x) the spherical Bessel function of order i. By substituting this result into equation (g) we can 
easily solve the Gaussian integration over the Wi. In order to perform the remaining integrations, it is convenient to 
introduce the new variables q and Q = qi + q 2 - In this way, after some algebra, we finally obtain: 

1 
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and £>(a;) represents the Dawson's integral ^ (see, e.g., Abramowitz & Stegun, 1968). It is straightforward to show 
that for t\ = i 2 the previous formula reduces to the usual expression for the mass two-point correlation function in 
ZA (Bond & Couchman 1988, Mann, Heavens & Peacock 1993, Schneider & Bartelmann 1995). 

We numerically evaluated the two-point correlation function f(r, t) = £(r,t,t) employing a COBE normalized 
standard CDM linear power spectrum (with density parameter SI = 1 and h — 0.5). We used the transfer function of 
Bardeen et al. (1986) while the normalization to the four-year COBE DMR data is given in Bunn & White (1997) and 
corresponds to ag = 1.22. As already noted by MHP, the small scale behaviour of the resulting correlation function 
depends on the value assigned to the truncation radius, Rf, defined in equation (^) (see Fig. |l|). If Rf is very small, 
then shell crossing will not be suppressed and £(r) will show an unusually flat behaviour. On the contrary, if Rf is 
too large, the smoothing procedure will remove an important contribution to the power spectrum, causing again too 



t It is worth stressing that when ip± < ip« , in order to avoid a complex argument for the Dawson's integral, it is convenient 
to express the integrand in equation Jll| ) in terms of exponentials and error functions (see also the discussion in Schneider & 
Bartelmann, 1995). However, since for the CDM spectrum (the only one considered in our analysis) tp^ is never larger than ipj_, 
we preferred to write the solution using D(x). 
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Figure 2. Redshift evolution of the mass two— point correlation function obtained using ZA to evolve a linear CDM spectrum. 

low a correlation. Therefore we need a criterion to select Rf. Since our main purpose is to compare the clustering 
amplitudes predicted by ZA with those extracted from the scaling ansatz of JMW, we can choose Rf so as to optimize 
the agreement between the respective correlation functions. Anyway, we find that this method conforms quite well 
to a simpler one already used by MHP: the best Rf is the one that maximizes £(r,Rf) on small scales. Strictly 
speaking, the optimal smoothing radius depends on the scale selected for maximizing the correlation: the smaller 
is r the larger comes out Rf (we find that the difference between the smoothing lengths obtained by maximizing £ 
at r — 0.1/i -1 Mpc and at r = l/i _1 Mpc roughly amounts to 0.2/i _1 Mpc and remains nearly constant by varying 
erg). However, the effect of this discrepancy on the correlation evaluated on larger scales is indeed minimal. Following 
Schneider & Bartelmann (1995), we select r — l/i -1 Mpc as the scale at which we require ((Rf) to be maximal. As 
previously stated, the optimum filtering length increases as the field evolves; the dependence of the best Rf on erg is 
almost linear and for og > 0.3 (that in our model corresponds to z ~ 3) it can be approximated by: 

Rf(a 8 ) = (3.16 ct 8 -0.65) /i _1 Mpc . (13) 

The redshift evolution of the correlation function is shown in Fig. ^. As expected, on scales that are not affected 
by shell crossing (r > Rf), £(r,z) steepens with decreasing z. Moreover, we note that the first zero crossing radius 
of £(r, z) increases as time goes on (see also Bond & Couchmann 1988). A similar pattern has been noticed by 
Coles (1990) and by Baugh & Efstathiou (1994) in the context of second-order Eulerian perturbation theory. The 
displacement of the first zero crossing of £ as a function of time is plotted in Fig. ^. Measuring the degree of dynamical 
evolution of the density field through erg, this shifting can be described with good approximation by the function: 

rocM - rll ~ 5.3a^ 5+0 - 1/<T8) fc~ l Mpc (14) 

where we denoted by roc the scale at which the correlation function crosses for the first time the zero-level and by 
t 1 qq its linear counterpart. It would be interesting to compare this result with the predictions of second-order Eulerian 
(and Lagrangian) perturbation theory and of other dynamical approximations. 



4 COMPARISON WITH THE SCALING HYPOTHESIS 

The analysis of a large set of numerical simulations suggests that, in hierarchical models, the non-linear two-point 
correlation function, £(r, z), can be related to the linear one, £h(r, z), through a simple scaling relation (Hamilton 
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Figure 3. The first zero crossing radius of the correlation function is plotted against 1/(78 • The circles represent the results 
obtained using ZA, the dashed line is the fitting function given in the text while the dotted line shows the prediction of Eulerian 
linear theory. 



et al. 1991, Peacock & Dodds 1994, Jain, Mo & White 1995, Peacock & Dodds 1996). The main idea is that the 
action of gravity can be represented as a continuous change of scale or, better, that the 'flow of information' about 
clustering propagates along the curves of equation: 



where £(r, z) represents the average correlation function within a sphere of radius r 



(15) 



(16) 

and ro is a sort of Lagrangian coordinate determining a 'conserved pair surface' (Hamilton et al. 1991, Nityananda 
& Padmanabhan 1994). In fact, by definition, the average number of neighbours of a particle contained within a 
spherical volume of radius r at the linear stage (when ( < 1) equals the average number of neighbours inside a 
sphere of radius r in the evolved field. 

Here, we want to compare the results obtained in the previous section, using ZA, with the predictions of the 
scaling ansatz (hereafter SA) formulated in the version of JMW: 



£,{r,z)=B(n c «)F 
with 
F(x) = 



B(n eS ) 



x + 0.45a; 2 - 0.02a 6 + 0.05a; 6 



1 + 0.02a; 3 + 0.003a; 9 / 2 



B(n cfi ) 



^ 3 + n eff y 



n e ft(z) = 



dlnP(fc) 



din A; 



(17) 



(18) 



i=nl( z ) 



where denotes the radius of the top-hat window function in which the rins linear mass fluctuation is unity. 
However, it would be useless to perform the comparison between the spherically averaged correlation functions since, 
on small scales, £ (r) obtained using ZA is seriously affected by shell crossing and the computation of £ requires an 
integration starting from r = 0. For this reason we prefer to use directly £(r). The two-point correlation function 
deriving from the ansatz of JMW can be obtained performing a simple differentiation: 



t(r,z) = 



[l + B(n cS )F(X)]F'(X)A^(r (h z) 
1 + B(n ea )F(X) - F'(X)A^ L (r , z] 

with F'{x) = dF/dx and 



+ B(n eB )F(X) 



X = 



$L(r ,z) 
B(n cB ) 



(19) 
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Figure 4. Comparison between the mass autocorrelation functions computed for a CDM model by using: the Zel'dovich 
approximation (ZA), the scaling ansatz of Jain, Mo & White (JMW) and Eulerian linear theory (ELT). The linear power 
spectrum extrapolated to the present epoch (z = 0) is normalized to match the COBE DMR data (trg = 1.22). 

Afr(ro,z)=&(r ,z)-fL(ro,z) = ^y J k 2 P(k)[j (kr ) - —jt(kr )]dk . (20) 

We evaluated the correlation function given in equation ( |l9| ) using a COBE normalized, linear CDM spectrum. In 
Fig. |^ we plot the result obtained at z — with the corresponding one achieved by using ZA. For comparison we also 
show the prediction of Eulerian linear theory. The agreement between ZA and SA is remarkable on mildly non-linear 
scales (4/i _1 Mpc^ r'tS 20/i _1 Mpc) and on completely linear scales (r > 50 /i _1 Mpc). For example, at r = 5/i _1 Mpc, 
linear theory overestimates the correlation of JMW by 82%, ZA underestimates it by 2% while the accuracy of the 
JMW fit is about 15 — 20%. However, we find that in the interval 20 /i _1 Mpc^ rv$ 50 ft _1 Mpc ZA predicts more non- 
linear evolution than SA (for example the roc obtained by using ZA is larger than the one determined through SA). 
In order to consider a less evolved field, in Fig. [| we repeat the comparison using the correlation functions evaluated 
at z = 1. Now, the main item to note is that the JMW result matches the linear solution on scales (r ~ 10/i _1 Mpc) 
that, according to ZA, are already involved in non-linear phenomena. 

In any case, we do not know the accuracy of the scaling hypothesis on large scales. In fact, the function F(x) 
is obtained by requiring the resulting £(r) to reproduce the linear behaviour where £l —*■ and, simultaneously, to 
approximate properly the correlation function extracted from N-body simulations. However, in order to achieve a 
detailed description of non-linear scales, JMW used a relatively small box to perform their simulations. Therefore, 
imposing the match to linear theory on large scales, without having any constraint from numerical data on quasi- 
linear scales, could seriously alter the accuracy of F(x). This probably implies that the JMW fitting function could 
be improved on large scales. Our conclusion is shared by Baugh & Gaztanaga (1996, hereafter BG), who tested the 
scaling ansatz for the evolution of the power spectrum against the results of 5 N-body simulations performed within 
a 378/i _1 Mpc box. Indeed, they found that the JMW formula gives a relatively poor description of the large-scale 
behaviour even though the agreement between the spectra remains always within the quoted 20 % accuracy. 

By using the output of their simulations, BG proposed a new scaling formula calibrated on large scales. As 
initially suggested by Peacock & Dodds (1994), the analytic expression of this SA concerns the dimensionless power 
spectrum, A 2 (fc, z) — k 3 P(k, z)/2n 2 (i.e. the contribution to the the variance of the density contrast per bin of Ink), 
while, following JMW, it takes account of a spectral dependence of the transformation: 

fc L = [l + A 2 (k,z)]~ 1/3 k (21) 




A 2 (fc, 



0(n e ff) 
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Figure 5. As in Fig. 4, but at z = 1 (erg = 0.61). 



where 



1 + 0.598s - 2.39a; 2 + 8.36x 3 - 9.01t 3 ' 5 + 2.895x 4 



1/2 



(3{n e 



1.16 



^ 3 + n cff \ 



1/2 



(22) 



1 - 0.424s + [2. 895/(11. 68) 2 ]a; 3 / ' <~ — y 3 

and the subscript L marks linear quantities. The function f(x) has been obtained by matching the power spectrum 
in the simulations at erg = 1, with an accuracy of 5%, over the range 0.02 /iMpc -1 < k < 1.0/iMpc" 1 and by forcing 
the fit to have the asymptotic form f(x) — > 11.68 x 3 ^ 2 when x — > oo (Hamilton a/. 1991). The two-point correlation 
function is related to A 2 (k,z) through the Fourier relation: 

dk 



£(»-,*) 



A'(M)i„(*r) 



(23) 



In Fig. |^ we compare the correlations obtained by using ZA and the JMW formula with the results of the scaling 
ansatz by BG: we are considering a standard CDM linear spectrum at the epoch in which as = 1.22. We immediately 
note that using larger simulation boxes to calibrate the SA allows a better determination of the correlation function 
for r> 20 h^Mpc. In fact, we find that the correlations obtained with ZA and with the BG formula agree by better 
than 20% for r > 4.6 h~ 1 Mpc (with the exception of a very small r-interval centred in the first zero crossing of £) while 
the discrepancy between ZA and the JMW ansatz is less than 20% over the ranges 4.1h~ 1 Mpc < r < 18.3/i _1 Mpc 
and r > 49.6/i _1 Mpc. Similar patterns are obtained considering different values of as- This shows that the BG fit, 
that has been calibrated against large box CDM simulations, gives also a very good description of the mass clustering 
predicted by ZA on intermediate scales. In any case, as expected, the JMW formula is sensibly more accurate for 
5/i -1 MpcJ; 15/i _1 Mpc where the BG predictions grow worse as as assumes values significantly larger than 1. 

On the other hand, it would be interesting to check the reliability of ZA and second-order Eulerian perturbation 
theory by directly comparing their predictions on these scales. Bond & Couchmann (1988), studying the weakly 
non-linear evolution of the CDM power spectrum, found remarkable agreement between the two approximations. 
Moreover, Baugh & Efstathiou (1994) showed that second-order Eulerian perturbation theory can reproduce, at least 
qualitatively, the evolution of the power spectrum predicted by numerical simulations. However, Jain & Bertschinger 
(1994) found that the agreement between perturbation theory and N-body outcomes gets worse as the density field 
evolves. Besides, their results are inconsistent with the low-fc behaviour of the second-order Eulerian correction to 
the CDM power spectrum computed by Bond & Couchmann (1988), raising again the issue about the compatibility 
between ZA and perturbation theory. In a recent work concerning the evolution of scale invariant spectra, Scoccimarro 
& Friemann (1996b) showed that, if the spectral index n satisfies —3 < n < —1, Eulerian perturbation theory is 
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Figure 6. Mass two— point correlation functions at the epoch in which <rg = 1 obtained from a linear CDM spectrum evolved 
through the ZeFdovich approximation (ZA) and through the scaling ansatze by Jain, Mo & White (JMW) and by Baugh & 
Gaztanaga (BG). 

able to reproduce fairly well the power spectrum obtained though the scaling ansatz, while the one-loop perturbative 
version of ZA gives worse results. Anyway, Bharadwaj (1996a, b) pointed out that the effects of multistreaming on 
the correlation function cannot be studied perturbatively. This fact implies that our result, obtained considering the 
full Zel'dovich approximation, should be more reliable than any other achieved by adopting a perturbative version 
of ZA. In any case, it would be interesting to clarify to which extent ZA and Eulerian perturbation theory agree on 
large scales. 




5 THE CORRELATION OF HIGH REDSHIFT OBJECTS 

In this section, we study the evolution of the cross correlation function of the mass density contrast evaluated at two 
different times as defined in equation (|ll|). This quantity could play an important role in comparing the clustering 
properties extracted from deep redshift surveys to the predictions of theoretical models for structure formation. In 
practice, one always collects data on correlations in a finite redshift strip of his past light cone while the quantity £ (r, t), 
normally used in theoretical works, refers to objects selected on an hypersurface of constant cosmic time. Therefore, as 
far as one is considering a deep sample of cosmic objects, it is not correct to relate the observed clustering properties 
to f(r,t). This issue is addressed in detail by Matarrese et al. (1997, hereafter MCLM) who build a theoretical 
quantity that allows a direct comparison of model predictions to the observed correlations. Their approach can be 
divided into three steps: first of all they compute the redshift evolution of mass correlations, then they relate the 
clustering properties of cosmic objects to the matter distribution by means of a linear bias relationship and finally 
they convolve the result with the observed redshift distribution of the class of objects under analysis. By assuming 
that the effects of redshift distortions and of the magnification bias due to weak gravitational lensing are negligible 
and by considering isotropic selection functions, MCLM showed that the theoretical estimate for the observed two- 
point correlation function can be formally expressed as an integral over z\ and 22 of the function f (r, Z\,Z2) weighted 
by geometrical factors and effective bias parameters (all dependent on z\ and 22). Different classes of objects are 
selected by changing the amplitude and the redshift dependence of the effective bias. However, in the absence of a 
model for the evolution of the cross-correlation, only assuming that the above mentioned integral is dominated by 
the contribution of objects whose redshifts are nearly the same, can one estimate the observed correlation function 
deriving from a particular scenario of structure formation. In this way, one is allowed to replace £(r, z\, 22) with £(r, 2), 
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Figure 7. Cross correlation between the density contrast field evaluated at two different redshifts vs. comoving separation. 

where z is a suitably defined average between Zj and 22 that, for simplicity, MCLM identify with z = (z\ -\-z%)j1. This 
is a crucial approximation, as it allows MCLM to use the JMW ansatz to compute the non linear mass correlation 
function (there is no known scaling ansatz for f(r,2i, 22)). However, as shown in the previous paragraphs, ZA allows 
the computation of £(r, 21,22) so that we are able to compute the theoretical estimate for the observed correlation 
function by using both the complete and the approximated formulae given by MCLM (respectively their equations 
15 and 18). Therefore we can check here, within the validity of ZA, the reliability of the approximation introduced by 
MCLM. Large discrepancies between the exact and the approximated correlations would obviously invalidate their 
whole analysis and consequently also their complete formula for £ D b s would be unutilizable. On the other hand, if 
the approximated correlation function turns out to reproduce accurately the complete one, MCLM formulae could 
represent an important tool to disprove cosmological models in the light of present and future observations. 

In order to compute £(r, 21,02) using equation (|ll|), we truncated the linearly extrapolated power spectrum 
b(z\)b(z2)P(k) according to the prescription: 



where Rf(z) represents the optimum filtering length for the density field at redshift z, determined by following the 
method described in Section 3. On small scales, the correlation functions that we obtain opting for this truncation 
procedure appear much more flattened than those computed at a single time. The evolution of £ (r, 21, 22) as Z2 changes 
is shown, for a CDM model, in Fig. [?]. It is evident that even though the correlation decreases as 22 grows, its decay 
is very slow. Actually, the ratios between the correlations computed at the same r, for different pairs of redshifts, are 
very similar to the predictions of linear theory. We find that the redshift evolution of the cross correlation function 
can be approximately described by the relation: 



where the quantity s = r/roc{z) is introduced in order to take into account the shifting of the first zero crossing of 
£(r, z) and O(x) is the Heaviside step function. Moreover, the first zero crossing radius of £(r, zi, Z2) is nearly given by 
the geometric average of roc(2i) and roc(z2). For s > 0.1 equation (|2^), which is meaningful up to the scale at which 
the first of the two £(s, z) reaches its second zero crossing, reproduces £(s, Zi, 22) with an accuracy of ~ 5%. Anyway, 
for 2, the usual relation £(r, 21,22) — [£(r, 2i)£(r, 22)] 1//2 sign[£(r, 21)] deriving from linear theory is preferable. 
We can now check the accuracy of the approximation introduced by MCLM that consists in computing the 




(24) 



21, 22) [£(s, 2 2 )] 1/2 [1 - 20(s - 1)] 



(25) 
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Figure 8. The observed two-point correlation function computed using equation (26) for a CDM model with [2 m i„,x ma x] = 
[0, 2]. For comparison, the corresponding f (r, z) evaluated for z = 0, 1, 2 are plotted. 

theoretical estimate for the observed correlation function by replacing zi, 22) with £(r, z), where z = (zi + £2)/2, 
in the appropriate formula. For simplicity (and in order to isolate the phenomenon of clustering evolution) we will 
assume no bias, no selection effects and a constant comoving number density in an Einstein-de Sitter universe. In 
this case, equation 15 of MCLM reduces to: 



where we denoted by ^obsO", 2 mm, Zmax) the (ensemble averaged) theoretical estimate for the two-point correlation 
function measured by an observer that acquires data from the region of his past light cone corresponding to the 
redshift interval [z m in, Zmax]- 

Considering only the linear evolution of density fluctuations, £(r, zx, 22) = £(r, 0, 0)/[(l+2i)(l+22)], the integrals 
contained in equation (^) can be analytically performed. In this case, the quantity £ a be(r,zx, Z2)/£(r, 0, 0) does not 
depend on r; for example we obtain £obs(r, 0, 2)/£(r, 0, 0) ~ 0.224 and £obs(r, 0, l)/£(r, 0,0) ~ 0.375. In this regime, 
we find that the approximation for £ b s introduced by MCLM is accurate to 2 — 3%. 

In order to extend our analysis also to the mildly non-linear evolution, we numerically computed £ b s by using 
the cross correlation given in equation (|ll|). The result obtained for [2 m in,Zmax] = [0,2] is shown in Fig. []: also in 
this case £ b s looks like the usual correlation function evaluated at some intermediate redshift. We then tested the 
accuracy of the above mentioned simplified scheme for the computation of £ b s , finding good agreement between the 
exact and the rough estimates (excluding a small neighbourhood of the zero-crossing radius of £(r, 21,02), where 
the approximated method breaks down, we find a maximum discrepancy of 6% for [z m i n , z max ] = [0,2] and of 3% 
for [zmin^max] = [0,1]). Anyway, the simplified procedure to compute £ bs can be further improved: adopting a 
different way of performing the average between redshifts, namely 1 + z = [(1 + 2i)(l + Z2)] 1 ^ 2 , ensures more accurate 
predictions (in this case the maximum error is always of the order of 1%). Probably this higher precision is due to 
the fact that we are considering mildly non-linear scales and the latter approximation gives exact results for linear 
evolution. 




(26) 
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6 SUMMARY 

In this paper, we have studied in detail the evolution of the mass two-point correlation function by describing the 
growth of density perturbations through ZA. Our motivations were originated by the well known ability of ZA to 
reproduce the weakly non-linear regime of gravitational dynamics. On scales that are not affected by shell-crossing, 
we found that the correlation function steepens as the clustering amplitude increases. Moreover, we showed that 
non-linear interactions are able to move the first zero crossing of £ (r) and we gave a quantitative description of this 
shifting for a CDM linear spectrum. 

We then compared our results with the predictions of the scaling ansatz for clustering evolution formulated by 
JMW, obtaining remarkable agreement between the correlations on mildly non-linear scales and on completely linear 
scales. However, between these two regimes, the JMW prescription, which has been obtained requiring the resulting 
correlation to reproduce the linear behaviour on large scales, predicts smaller clustering amplitudes than ZA. We 
think that this disagreement is caused by the smallness of the box used by JMW to perform their N-body simulations. 
Actually, imposing to match the linear solution where £l — > 0, without having any constraint from numerical data on 
quasi-linear scales, could alter the accuracy of the fitting function that embodies the scaling ansatz. In connection with 
this hypothesis, we compared ZA predictions on correlations with the output of a different scaling ansatz calibrated 
against large box simulations by BG. In effect, on large scales, the BG formula agrees better with ZA, keeping the 
same accuracy of the JMW fit on intermediate scales. 

On the other hand, the reliability of ZA on these scales and for dynamically evolved fields (cts~ 1) should be 
verified by directly comparing its predictions with the results of other approximations and numerical simulations. 

Finally, we studied the evolution of the cross correlation between the density field evaluated at two different 
epochs and, adopting the method introduced by MCLM, we used our results to compute the theoretical prediction 
for the observed correlation function deriving from a deep catalogue of objects. In this context, we proposed a 
simplified procedure for the computation of £ obs that, at least for quasi-linear scales, significantly improves another 
approximation previously introduced by MCLM. This result confirms that the MCLM method can be used to make 
quantitative predictions about clustering evolution that find a direct observative counterpart in the analysis of deep 
surveys. 
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